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Abstract 

Along with tt° and n mesons, a resonance structure in the invariant mass 
spectrum of two photons at M 77 = 360 ± 7 ± 9 MeV is observed in the reaction 
dC — > 7 + 7 + X at momentum 2.75 GeV/c per nucleon. Estimates of its width 
and production cross section are T = 63.7 ± 17.8 MeV and ct 77 = 98 ± 24+^ ub, 
respectively. The collected statistics amount to 2339 ±340 events of 1.5- 10 6 trig- 
gered interactions of a total number ~ 10 12 of dC-interactions. This resonance 
structure is not observed in pC collisions at the beam momentum 5.5 GeV/c. 
Possible mechanisms of this ABC-like effect are discussed. 
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1 Introduction 

Dynamics of a near-threshold production of the lightest mesons and their interactions, 
especially the pion-pion interaction, are of lasting interest. A good understanding of the 
pion-pion scattering is essential as it provides a test for the chiral perturbation theory 
and the information about quark masses and the chiral condensate. The two-photon 
decay of light mesons represents an important source of information. In particular, the 
77 decay of light scalar mesons was considered as a possible tool to deduce their nature. 



Also, the scalar-isoscalar sector is under much debate presently since more states are 
known (including possible glueball candidates) than can be fitted into a single multiplet. 
Unfortunately, the existing experimental information from irir scattering has many 
conflicting data sets at intermediate energies and no data at all close to the threshold 
region of interest. For many years this fact has made it hard to obtain the conclusive 
results on im scattering at low energies or in the sigma region. 

The so-called "ABC effect" is among the oldest and still puzzling problems. Al- 
most fifty years ago Abashian, Booth, and Crowe [1] first observed an anomaly in the 
production of pion pairs in the reaction dp — > 3 He+27r = 3 He+A°. This anomaly or 
ABC effect stands for an unexpected enhancement in the spectrum of the invariant 
7T7T mass at masses of about 40 MeV higher than 2m w . The subsequent experiments 
dp ^ 3 Re+X 0, pn -»• d + 2tt SI E], dd ^ 4 He+A [7] , np -> d + 2vr with neutron 
beams [H [9] and even np —* d + 77 [10] independently confirmed this finding. This 
anomaly was also observed in the photoproduction of pion pairs, ^p — > p + X° [ITJ 12] • 
It was revealed that the ABC effect is of isoscalar nature since a similar effect was 
not observed in the pd — > H 3 + X + reaction. The peak positions and widths vary for 
different bombarding energies and observation angles. Initially, the low-mass enhance- 
ment was interpreted as being caused by an unusually strong s-wave nir interaction 
or as an evidence for the a meson existence [TJ which shortly before was suggested 
by Johnson and Teller to provide saturation and binding in nuclei \13\. It is usually 
accepted now that this enhancement is not an intrinsic two-pion property since there 
is no resonance structure in the ttit scattering amplitude in this energy range. So any 
interpretation of the ABC as a real resonance is very much in doubt (for example, see 
discussion in [8]). It is generally believed that a system like that has to be associated 
with two nucleons when two pions (both must be present) are rescattered off them or 
both nucleons participate in elementary pp — > it + X reactions (predominantly via A 
formation). Actually, the origin of the ABC effect must be looked for in the formation 
of light nuclei at intermediate energies (for a review see Ref . [H] ) . 

The presented complicated situation is reflected in the Particle Data Group (PDG) 
table [15] where the values quoted for the sigma mass and width, based on both the 
pole position and the Breit-Wigner parameter determinations, are very widely spread. 
The estimated mass and half-width are 



However, during the last years, the chiral perturbation theory and Roy equations led to 
an accurate description of ttit scattering at low energies and the precise determination 
of the mass and width of the a resonance [16]: 



All experiments conducted on the ABC issue with the exception of low-statistics 
bubble- chamber measurements [2J Ej were inclusive measurements carried out prefer- 
entially with a single-armed magnetic spectrograph for detection of the fused nuclei. 
They allow one to find the two-pion invariant mass through the missing mass. Very 
recently, exclusive measurements of reactions pd — > p + d + 7r° + 7r° and pd — > 3 He+7r + 71 
have been carried out with complete kinematics in the energy range of the ABC effect 
at CELSIUS using the An WAS A detector [13 [T5]. The importance of the strongly 
attractive AA channel was noted. Surprisingly, the basic pp —> pp7r°7r° reaction in 



m a - TV/2 = (400 - 1200) - z(250 - 500) MeV. 



(1) 



m a = Ultl 6 MeV; Y a /2 = 272± 9 12 5 MeV. 



(2) 



the AA region also shows an ABC-like 1ow-7T7t mass enhancement which deserves spe- 
cial attention. It confirms the earlier result in [19] where the analyzing powers and 
cross sections for the ABC enhancement production were measured for the reaction 
pp — * p + p + X° in the missing mass range < M x o < m v . However, the big differ- 
ence was observed in the width of the resonance cross section and it was concluded that 
the observed width in the isoscalar channel is not obviously just a simple result of the 
binding between the two A states. It rather signals more complicated configurations 
in the wave function of the intermediate state, as would be expected for a nontrivial 
dibaryon state [IB] . 

This work aims to study whether this 1ow-7T7t mass anomaly can survive in heavier 
systems in the 77 channel. The paper is organized as follows. After a brief description 
of the experiment and experimental setup, the structure of measured invariant mass 
spectra of photon pairs is analyzed in Sect. 2. As a cross-check, a similar analysis is 
carried out in Sect. 3 but within the wavelet method. To elucidate the nature of the 
peak at M 77 pa (2 — 3)m n , different mechanisms of the observed 77 pair enhancement 
are discussed in Sect. 4. In Sect. 5, experimental estimates for production cross sections 
and widths of rj mesons and hypothetical R resonance are given. The main inferences 
of the paper are presented in the Concluding Section. 

2 Experiment 
2.1 General layout 

The data acquisition of production of neutral mesons and 7-quanta in pC and dC 
interactions has been carried out with internal beams of the JINR Nuclotron [201 EI] • 
The experiments were conducted with internal proton beams at momentum 5.5 GeV/c 
incident on a carbon target and with 2 H, 4 He beams and internal C-, A1-, Cu-, W- 
, Au-targets at moments from 1.7 to 3.8 GeV/c per nucleon. For the first analysis 
the data with the maximal statistics, pC- and dC-interactions, were selected. The 
first preliminary results on <iC(2 AGeV) collisions were reported in [21] where some 
indication on unusual structure in the photon-photon invariant spectrum has been 
obtained. 

The presented data concern reactions induced by deuterons with a momentum 2.75 
GeV/c per nucleon and by protons with 5.5 GeV/c. Typical deuteron and proton 
fluxes were about 10 9 and 2- 10 8 per pulse, respectively. The electromagnetic lead glass 
calorimeter PHOTON-2 was used to measure both the energies and emission angles of 
photons. The results obtained in earlier experiments with this setup are published in 
|22j . The experimental instrumentation is schematically represented in Fig{TJ 

The PHOTON-2 setup includes 32 7-spectrometers of lead glass and scintillation 
counters S\ and S2 of 2 x 15 x 15 cm 3 used in front of the lead glass for the charged 
particle detection [22| 1231 124] . 

The center of the front surfaces of the lead glass hodoscopes is located at 300 cm 
from the target and at angles 25.6° and 28.5° with respect to the beam direction. This 
gives a solid angle of 0.094 sr (0.047 sr for each arm). Some details of the construction 
and performance of the lead glass hodoscope are given in TablJU The internal target 
consists of carbon wires with the diameter of 8 micron mounted in a rotatable frame. 
The overall construction is located in the vacuum shell of the accelerator. 
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Figure 1: The schematic drawing of the experimental PHOTON-2 setup. The Si and 
5*2 are scintillation counters. 

Before the experiment the energy calibration of the lead glass counters has been 
carried out with 1.5 GeV/c per nucleon deuteron-beam at the JINR synchrophasotron 
|23j . The long-term gain stability was continuously monitored for each of the lead glass 
modules with 32 Nal(Tl) crystals doped with 241 Am sources. 

The modules of the 7-spectrometer are assembled into two arms of 16 units. These 
modules in each arm are divided into two groups of 8 units. The output signals of 
each group from 8 counters are summed up linearly and sent to the inputs of four 
discriminators (-Dj). In this experiment all the discriminator thresholds were at the 
level of 0.4 GeV. Triggering takes place when there is a coincidence of signals from 
two or more groups from different arms: (Di + D2) x (D3 + D4). In realizing the 
trigger conditions the amplitudes of all 32 modules were recorded on a disc. The dead 
time of data acquisition is about 150 fis per trigger. The mean rate of triggering was 
about 330 and 800 events per spill in dC and pC reactions, respectively. Duration of 
the irradiation cycle is 1 second. Totally about 1.52 x 10 6 and 1.06 • 10 6 triggers were 
recorded during these experiments. 

2.2 Event selection 

Photons in the detector are recognized as isolated and confined clusters (an area of ad- 
jacent modules with a signal above the threshold) in the electromagnetic calorimeter. 
The photon energy is calculated from the energy of the cluster. Cluster characteristics 
were tested by comparison with Monte-Carlo simulations of electron-photon showers 
in Cherenkov counters by means of the program package EMCASR [25] . The results 
obtained earlier with extracted ion beams of the JINR Synchrophasotron have demon- 
strated a good agreement between experimental and simulated data [22]. Assuming 
that the photon originates from the target, its direction is determined from the geomet- 



Table 1: The basic parameters of the lead glass hodoscope 



Number of lead glasses 


oZ lr-1, total wignt 14zz Kg 


Module cross section 


r=9 cm of insert circumference 


Module length 


oo cm, 14 rt. Li. 


T~\ Q T 1 Q 1 T*OC/~\ 1 11 "t" 1 /~\T1 

DJJclllcll lebUlULlUIl 


3.2 cm 


Angular resolution 


0.6° 


Energy resolution 


(3.9/ VE + 0.4)%, E [GeV] 


Gain stability 


(1-2)% 


Dynamic range 


50 MeV - 6 GeV 


Minimum ionizing signal 


382 ±4 MeV of the photon 




equivalent 


Total area 


0.848 m 2 



rical positions of constituent modules weighted with the corresponding energy deposit 
in activated modules. 

The invariant mass distributions of photon pairs (from different arms of the spec- 
trometer) are shown in Figj2j The dominant part of distributions (two upper panels) 
comes from the 7r° — > 77 decay. These photons in combination with others result in a 
huge background which mask the expected rj — > 77 decay. Other sources of background 
are charged particles as well as neutrons and particles from a general background in 
the accelerator hall. Contributions of the given sources were estimated by special 
measurements with and without veto-detectors SI and 5*2 and by comparison of data 
obtained at different beam intensities. The contribution to the total combinatorial 
spectrum of charged particles, neutrons and particles from a general background in the 
accelerator hall is less than 10% and becomes negligible (< 1%) after subtraction of 
mixing event background (see below). As follows from FigJ2l the high-energy cut of 
photons, Ery > 100 MeV, allows one to improve the signal-to-background ratio. The 
contribution of the general background in the experimental hall was estimated from 
the measurements with the empty target, see Figj3j Two runs carried out for 125 
accelerator cycles of the deuteron beam (about 10 11 deuterons in every run) result in 
Nyy = 117428 and 338 photon pairs for the case with and without target, respectively. 
So this source contributes less than 1% and is quite smoothly distributed with respect 
to M 77 . 

To see a possible structure of the invariant mass spectra, a background should be 
subtracted. The so-called event mixing method was used to estimate the combinatorial 
background: a photon in one event from the first arm is combined with a photon in 
other events from the second arm. This background was subtracted from the invariant 
mass distributions (see bottom panels in FigfS]). The background normalization was 
carried out in two steps. First, the background is normalized to the total pair number 
(see [21]). Naturally that at the event mixing the resonance maxima are not repro- 
duced disturbing the overall normalization. At the second step this shortcoming of the 
background estimate is corrected by an auxiliary factor K norm (which differs by few % 
from 1) obtained by iterating treatment of the resonance contribution to the spectrum. 
The negative values in the high-mass range arising at the subtraction of the mixing- 
event background are caused by the energy conservation law which may be violated 
when 7-quanta are taken from different mixing events. To decrease its influence the 
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Figure 2: (color online) Invariant mass distribution of 77 pairs from the reaction 
dC — > 7 + 7 + X at 2.75 GeV/c per nucleon for two values of the cut energy of photons. 
The top shaded histograms show the background contribution. The bottom histograms 
are invariant spectra after the background subtraction. The auxiliary normalization 
factor K norm (see below) is 0.9947 for left figure (the cut energy is 50 MeV) and 



0.993 for right figure (the cut energy is 100 MeV). 



energy sum of 7-quanta in both individual events and mixing 77 pairs is restricted (see 
below criterion (3)). A clear peak from the rj decay and a remnant of the suppressed 
7T° resonance are clearly observable. Note that between them there is some additional 
structure which will be clarified below. 

Systematic errors may be due to uncertainty in measurements of 7 energies and 
inaccuracy in estimates of the combinatorial background. The method of energy re- 
construction of events is described in detail in Refs. [22J Possible overlapping 
effect was studied at higher intensities in CC collisions at the 3.7 GeV per nucleon 
beam energy. It may result in about 6% displacement of the total mass spectrum. 
General influence of the overlapping effect in this experiment is very small. One of the 
criteria of accuracy of energy reconstruction is the conformity of the peak positions 
corresponding to the known particle mass values. As is seen in Fig{21 the position of 
peaks corresponding to rj- and 7r°-mesons is in reasonable agreement with the table 
values of their masses. A more precise determination of the position of peaks requires 
minimization of systematic errors in describing the background which arise, in par- 
ticular, due to the violation of the energy-momentum conservation laws in selecting 
7-quanta by random sampling from different events (see also below). 

The selection criteria can be made harder by imposing additional trigger conditions. 
For a background suppression and minimization of systematical errors due to violation 
of conservation laws the following selection criteria were used: 
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Figure 3: (color online) Invariant mass distributions of 77 pairs in two different runs of 
measurement under condition E 1 > 50 MeV: with the empty target (dashed histogram) 
and with the internal carbon target (solid histogram) in the reaction dC = 7 + 7 + X 
at 2.75 GeV/c per nucleon. 



(1) the number of photons in an event, N 1 = 2; 

(2) the energies of photons, E 1 > 100 MeV; 

(3) the summed energy in real and random events < 1.5 GeV. The criterion (1) 
suppresses the combinatorial background and minimizes systematical errors arisen due 
to the violation of photon topology at event mixing because events with a different 
number of photons iV 7 have different angular and energy distributions. As noted above, 
the criterion (2) improves the signal-to-background ratio. The criterion (3) also allows 
one to minimize systematical errors at event mixing at an insignificant (~ 3%) loss of 
events. The result of this triggering is shown in FigJU Under the selected condition 
the 7T° peak is practically absent. Therefore, 7r°-mesons were mainly detected in events 
with N y > 2 (a minimal opening angle of the 7 pair detected by the setup equals 42°). 
In contrast, the r\ is seen very distinctly with the width defined by the experimental 
resolution in the mass. In addition, in this reaction dC — > 7 + 7 + X a pronounced peak 
is observed in the interval 300-420 MeV of the invariant mass of two-photon spectrum 
which will be named below the R-resonance. 

However, under similar conditions only 77 is seen in the pC collisions. The observed 
peaks were approximated independently by the Gaussian: 
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(3) 



The additional shift-parameter y is introduced in eq.([n]). The values of the obtained 
fitting parameters are given in TablJJl 

The signal-to-background ratios for the invariant mass intervals of 300-420 MeV 
and 480-600 MeV (the vicinity of the ?7-meson mass) are 2.5 ■ 10~ 2 and 1.4 ■ 10 _1 , 
respectively. For comparison, analogous values without the background suppression 
(without the selection criteria (l)-(3)) are (4.0 ± 1.4) ■ 10" 3 and 3.2 ■ 10" 2 . 
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Figure 4: (color online) Invariant mass distributions of 77 pairs satisfying criteria 
(1) — (3) without (upper panels) and with (bottom panels) the background subtraction 
obtained for the reaction dC — > 7 + 7 + X at 2.75 GeV/c per nucleon (left) and pC 
collision at 5.5 GeV/c (right), respectively. The curves are the Gaussian approximation 
of experimental points (see the text). The values of K norm are 0.958 for dC and 0.952 
for pC. 



Table 2: Fit parameters of the Gaussian distribution 





dC 


dC 




165 < M 77 < 435 MeV 


465 < M 77 < 825 MeV 


yo 


-1.94 ± 1.31 


-0.004 ± 0.046 


No 


2623 ± 472 


7329 ± 295 


wWsur, MeV 


41.3 ± 7.2 


38.5 ± 1.7 


M Q , MeV/c 2 


362.0 ± 6.9 


535.7 ± 1.9 


X 2 / degrees of freedom 


9.06/6 


7.38/9 




pC 


pC 




345 < M 77 < 645 MeV 


405 < M 77 < 645 MeV 


Vo 


1.44 ± 1.86 


1.30 ± 3.07 


No 


5283 ± 560 


4804 ± 673 


Wmeasur, MeV 


51.6 ± 4.1 


41.9 ± 4.2 


M , MeV/c 2 


541.5 ± 2.5 


536.6 ± 2.6 


X 2 /degrees of freedom 


16.10/7 


3.13/4 
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Figure 5: (color online) The invariant mass distributions of two photons for the opening 
angles 0.55 < cos6 77 < 0.65 (left) and 0.65 < cos6 77 < 0.75 (right) under the 
selection criteria (1) -j- (2). The values of K norm are 0.98 (left) and 0.97 (right). 



Thus, as follows from TablfS], the position and width of the peak corresponding to rj- 
meson are in good agreement with values from the PDG table (systematic errors do not 
exceed 2%) and the spectrometer mass resolution. The total number of detected events 
in the dC reaction for the 77-meson region 450-660 MeV after background subtraction 
is 7336 ± 284. 

To elucidate the nature of the detected enhancement, we investigate the dependence 
of its position and width on the opening angle of two photons and on their energy 
selection level. The results demonstrated in Figs. [5] and [6] show that both maxima 
survive and are located practically at the same values of M 77 . Pair distributions 
over the opening angle 77 for two intervals of the sum of two-photon energy 0.8 < 
E Xl + E 2l < 1.1 (right) and 1.1 < E\ n + E 2l < 1.5 (left) are displayed in Fig. 
For lower values of the sum energy the two peaks are seen where the peak at smaller 
opening angles corresponds to the R resonance while 77 mesons are emitted at larger 
angles around 77 ~ 60°. Harder energetic selection (left panel) leaves only the 77 
meson contribution. 

As is seen, the result of the changing of observation conditions is quite robust: the 
position and width of the observed peak remain almost unchanged in different intervals 
of both energies and opening angles of 7-quanta, namely, the mean peak position in the 
invariant mass distribution varies under different conditions in the range 348-^-365 MeV. 
The total number of detected events in the region 270-450 MeV (a summed number of 
pairs in the histogram in FigH]) after the background subtraction is 2339 ± 340. 

The resonance-like structure observed in invariant mass di-photon distributions is 
not visible in energy photon spectra. As shown in FigJH] the 7 energy spectra near 
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Figure 6: (color online) The invariant mass spectra of 77 pairs for the energy selection 
E y > 400 MeV under the selection criteria (1) -j- (2) , K norm = 0.973. 

the R resonance are quite smooth, well reproduced by model simulations and have a 
very similar shape in both arms. This point testifies also that there is no instrumental 
anomaly causing the R peak. 

Because of smallness of the signal-to-background ratio in the R-resonance mass 
range, high statistics is needed for observation. The only most statistically meaningful 
measurement of the invariant mass spectra was made by the TAPS-collaboration [26] . 
The closest reaction where the 77 production was studied is CC interactions at the 
kinetic energy Tc =2 AGeV [26]. In this experiment 13290 ± 340 eta mesons were 
measured but a resonance structure in the range of M 77 = 300 -j- 400 MeV has not 
been recorded. If the R/rj ratio is assumed to be the same as in the case of dC collisions, 
one may expect about 1800- (e^/e^) of the detected R — > 77 events, where e# and e v are 
the detection and selection efficiency for R and 77, respectively. Taking into account the 
systematical uncertainties, this estimated value is in the limits of (500^-3000) • (e^/e^) 
of the detected R — > 77 events (see below, formula (fTUl) ). The total number of detected 
events in this range of two-photon invariant masses is about 6. ■ 10 5 [26], so roughly the 
resonance structure should be observed at the level of ~ 1800- (e^/e^) ±800. To resolve 
this structure general statistics should be increased by an order of magnitude, at least. 
Note that these two experimental setups cover different rapidity regions, which was not 
taken into account in our estimate. As compared to the TAPS, the PHOTON-2 setup 
has a smaller angular acceptance but a better signal-to-background ratio. Thus, there 
is no discrepancy between our result and observation of no resonance structure by the 
TAPS [26] in the reaction close in the energy and mass numbers. 

An indication to a possible resonance structure in pp-collisions at the energies 
1.36 and 1.2 GeV was obtained by the CELSIUS- WASA collaboration [52]. However 
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Figure 7: (color online) Distribution of the opening angle of 77 pairs in dC collisions 
for the two selections of (E^ + E 2l ). Other conditions are the same as in FigJH 
Experimental (circles) and simulated (triangles) results are normalized to the same 
number of photon pairs. The values of K norm are 0.974 (left) and 0.99 (right). 



statistics in the invariant mass range 250-350 MeV is low: only about 200 77 pairs 
without background subtraction was found. 



3 Wavelet analysis 

Here we shall try to identify essential structures in the measured M 77 spectra without 
the background subtraction. A conventional method for such analysis is a wavelet trans- 
formation which is known as an efficient multiscale technique to reduce the presence 
of statistical noise and then extract physical parameters from the obtained smoothed 
form [27J [2S]- The one- dimensional wavelet transform (WT) of a signal f(x) has a 
biparametric form. This allows WT to overcome the main shortcomings of the Fourier 
transform such as nonlocality, infinite support and necessity of a broad band of fre- 
quencies to decompose even a short signal. The wavelet transformation changes the 
decomposition basis into functions which are compact into a time/space and frequency 
domain. The WT with the wavelet function ip of the function f(x) is defined by 
convolution as 

W4a, b)f = -L = 7 -L^ (—) f(x) dx (4) 

VW-00 vi G i v a 1 

with the normalization constant 

Cw, = 2tt J duo < 00, 

J \u\ 
—00 

where ?p(ou) is the Fourier transform of the wavelet ip(x). A scale parameter a char- 
acterizes the dilatation and b is the translation in time or space. In this respect the 
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Figure 8: (color online) Photon energy spectra from 77 pairs in the invariant mass 
interval AM 77 = 0.32 -h 0.4 GeV. Experimental (circles) and Monte-Carlo simulation 
(triangles) points calculated with inclusion of the R resonance formation are given 
separately for every spectrometer arm. Distributions are normalized to the same total 
number of events. 




wavelet function ip(t) is a sort of a "window function" with a non-constant window 
width: high-frequency events are narrow (due to the factor 1/a), while low-frequency 
wavelets are broader. The inverse transform is given by the formula 

da db 

W^(a,b) — — . (5) 
cr 

The wavelet ip exists if < 00. It holds, in particular, when the first (n — 1) 
moments are equal to zero 

00 

\x\ m ip(x)dx = 0, < m < n. (6) 

Due to freedom in the choice of the wavelet function ip, many different wavelets were 
invented [2§J ED]- We consider here only continuous Wavelets with Vanishing Moments 
(WVM) (see it Appendix A). The WVM family is called so because condition (jH]) 
always holds for it. One of WVM families is a set of Gaussian wavelets (GW) which 
are normalized derivatives of the Gauss function 

g(x;Ax ) = Aexp (~ ( * "J?^ ) • (7) 

In some cases continuous wavelets are more suitable to evaluate peak parameters. 
One of these cases arises when a peak in question has a Gaussian shape (El|. This 
makes it possible to use very simple analytical expressions, eq. (fT9l) . in the continuous 
Gaussian wavelet transform for Gaussian peaks. It gives us a remarkable advantage to 
calculate the peak parameters directly in the wavelet domain instead of the time/space 
domain without using the inversion. Moreover, in real cases, when our signal shape is 
close to a Gaussian one and is considerably contaminated by an additive noise and, 
in addition, is distorted by binning to be input to computer, one can also use the 
remarkable robustness of Gaussian wavelet filtering, as proved in 



Wavelet ability to separate signal components with different frequencies and po- 
sitions has attracted many physicists to use both discrete and continuous wavelets 
[33l [M]. Usually the conventional filtering approach is applied: a signal transformed 
by a wavelet undergoes an appropriate thresholding and then is restored by the in- 
verse transform. The image of the wavelet spectrum is used to obtain rough parameter 
estimations of wanted peaks of invariant mass spectra, as in [35] where the Mexican 
hat wavelet was used. In our paper, we apply the family of GW to look for peaks in 
question having a Gaussian shape (J7|). 
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Figure 9: (color online) The invariant mass distribution of 77 pairs (top) and the 
biparametric distribution of the GW of the 8-th order (bottom) for pC interactions. 
These events are selected under the same conditions as in FigJH but the distribution is 
obtained with an additional condition for photon energies E^i/E^ > 0.8 and binning 
in 2 MeV. 
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Figure 10: (color online) The same as in Figl9]but for dC collisions. 

The main idea of our approach is to transform the signal f(x) to the space of the 
corresponding wavelet and look there for a local biparametric area surrounding the 
wavelet image of our peak in question, drawing down all other details of the signal im- 
age, concerning noise, binning effects, and background pedestal. A particular example 
of the (^2(0, b) transform and other details of this analysis are given in Appendix A. 

The initial signal f(x), i.e. the M 77 distribution including the background is pre- 
sented in the upper panels of Figsj9] and [101 Here an additional condition on photon 



energies E^i/E^ > 0.8 (where E y i and E l2 are the smallest and the largest energy 
in the given photon pair, respectively) was applied to data to improve the signal-to- 
background ratio and avoid double-humped structure which will require higher orders 
of the expansion. The distributions look quite smooth due to the choice of a smaller 
step in M 77 which is needed to provide more points for the wavelet analysis. 

The wavelet transformation results are presented in the lower panels of Figs. M and 
[TU1 for pC and dC, respectively. Assuming symmetric signal we use the WVM of the 
8-th order to separate noise and see a peak structure. Attempts to apply wavelets of 
lower than 8-th order give worse results, perhaps, because of rather ragged signals. 
The arrows show an approximate location of the identified peaks. Due to the trigger 
condition, the distribution maximum (a) of photons from the ir° decay is shifted to 
M 77 ~ 125 MeV from the expected 135 MeV. It is in agreement with the initial 
distribution shown in the same figures. A huge peak from the background (b) dominates 
and its shape is quite close to the Gaussian. A photon peak from the r\ decay, (d), is seen 
quite distinctly at the proper place but its intensity is suppressed due to an additional 
condition E^i/E^ > 0.8. In the domain of interest there is some enhancement near 
M 77 ~ 370 MeV marked in the figures by (c). In the biparametric representation 
this spot has not a circle shape because even in the case of a coarse binning it is not 
well approximated by the Gaussian (see TablfJ]). In this respect, the (c)-peak in pC 
seems to be more pronounced than the appropriate one in dC but it follows from lower 
statistics in the pC case where a separate point may be better approximated by the 
Gaussian. Note that statistics in these cases differs by the factor of more than 3. The 
WVM analysis still reveals one more weak (c)-peak at higher M 77 . This is not very 
surprising since in coarse binning (see Figf5]) they were blurred but they are seen at a 
more strict selection (cf. FigsJH -j- ED- 

Therefore, the presented results of the continuous wavelet analysis with vanishing 
moments confirm the finding of a peak at M 77 ~ (2 — 3)m n in the 77 invariant mass dis- 
tribution obtained within the standard method with the subtraction of the background 
from mixing events. 

4 Data simulation 
4.1 About the model 

To simulate pC and dC reactions we use a transport code. At high energies it is the 
Quark-Gluon String Model (QGSM) [36] and at the energy of a few GeV the string 
dynamics is reduced to the earlier developed Dubna Cascade Model (DCM) [37] with 
upgrade of elementary cross sections involved [38] . 

The DCM divides the collision into three stages, well separated in time. During 
the first initial stage an intranuclear cascade develops, primary particles can scatter 
and secondary particles can re-scatter several times prior to their absorption or escape 
from the nucleus. At the end of this step the coalescence model is used to localize 
d, t, 3 He, and 4 He particles from nucleons found inside spheres with well-defined radii 
in configuration space and momentum space. The emission of cascade particles de- 
termines a particle-hole configuration, i.e., Z,A, and excitation energy that is taken 
as the starting point for the second, pre-equilibrium stage of the reaction, described 
according to the cascade exciton model [39]. Some pre-equilibrium particles may be 



emitted and this leads to a lower excitation of the thermalized residual nuclei. In the 
third, final evaporation/fission stage of the reaction, the de-excitation of the residue 
is described with the evaporation model. The last two stages are important for trig- 
gering the events. All components contribute normally to the final spectra of particles 
and light fragments; low-energy evaporated photons are not included into subsequent 
analysis. For relativistic energies the cascade part of the DCM is replaced by the 
refined cascade model, which is a version of the quark-gluon string model (QGSM) 
developed in [10] and extended to intermediate energies in [¥!] . The description of 
the mean-field evolution is simplified in the DCM in the sense that the shape of the 
scalar nuclear potential, defined by the local Thomas-Fermi approximation, remains 
the same throughout the collision. Only the potential depth changes in time, according 
to the number of knocked-out nucleons. This frozen mean-field approximation allows us 
to take into account the nuclear binding energies and the Pauli exclusion principle, as 
well as to estimate the excitation energy of the residual nucleus by counting the excited 
particle-hole states. This approximation is usually considered to work particularly well 
for hadron-nucleus collisions. 

The following 7-decay channels are taken into account: the direct decays of 7r°, 77, r( 
hadrons into two 7's, u — > 7r°7, A — > Nj and the Dalitz decay of 77 — > 7r + 7r _ 7, 
77 — > 7e + + e~ and tt° — ► / ~fe + + e~, the 77' — > p° + 7, the E — > A + 7, the irN and 
AfTV-bremsstrahlung. One should note that in accordance with the recent HADES 
data [12], the pn-bremsstrahlung turned out to be higher by a factor of about 5 than a 
standard estimate and weakly depends on the energy. This finding, being in agreement 
with the recent result of Ref . [13] . allowed one to resolve the old DLS puzzle [HJ. This 
enhancement factor is included in our calculations. 




Figure 11: The proton energy dependence of the double differential cross section for 
the 77 production in pC collisions. Experimental points are from [15] . 




Figure 12: (color online) Transverse momentum distributions of 7r° and 77 in the middle 
rapidity range from CC collisions at different energies. Experimental points are from 
the TAPS Collaboration [26]. 

As a model test, in FigfTHthe excitation function for the 77 production is shown for 
the pC collisions. The model describes correctly a fast increase of the 77 yield near the 
threshold where the cross section is changed by two orders of magnitude. 

The transverse momentum distributions at the mid-rapidity are presented in FigJT2l 
for 7T° and 77 produced in CC collisions at three bombarding energies. The model results 
are in good agreement with the TAPS experiment (26] for both neutral pions and eta 
mesons. So this gives us some justification for application of our model to analyze 
neutral particle production in the reactions considered. 

4.2 Analysis of the obtained data 

The model described above is implemented for describing the measured distributions 
with careful simulations of experimental acceptance. The total statistics of simulated 
events amounts to about 10 9 here and in every case below. As is seen from FigJT3l the 
model reproduces quite accurately the observed 77 peak in the invariant mass distribu- 
tion of 7 pairs but there is no enhancement in the region of the i?-resonance. 

The mass of the R resonance is slightly above 2m n . As was noted, there is no 
particularities in the phase shifts for the irir scattering. We remind that according to 
the PDG table [15], the closest-in-mass hadron in this energy range is the /o(600)- 
meson (or the cx-meson) with the extended mass and very large width (see eq.([T])); 
however, the recent analysis gives a more defined a mass and a more narrow width (see 
eq.©)- 

To see whether such a resonance structure would be created due to a nuclear in- 
teraction and can survive in the strict experimental conditions, we artificially simulate 
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Figure 13: (color online)The invariant mass distributions of 77 pairs from the dC (left) 
and pC (right) reactions after background subtraction. Both experimental (circles) and 
simulated (triangles) points are obtained under the same PHOTON-2 conditions. 
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Figure 14: Invariant mass distribution of pion pairs from tttt interactions in pC (left) 
and dC (right) collisions. The contributions resulting in the /^-resonance and formation 
of p,uj,r]' mesons are shown separately. 
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Figure 15: (color online) Invariant mass distributions of 77 pairs from the dC and 
pC reactions after background subtraction. Both experimental (circles) and simulated 
(triangles) points are obtained under the same conditions. The contribution of photons 
from the R decay is shown by the solid line. 



production of the i?-resonance and follow its fate in the course of a nuclear collision. It 
is assumed that the putative R resonance can be created in every 7r + 7r~ or 7r°7r° inter- 
action if its invariant mass satisfies the Gaussian distribution with the observed 
parameters (see Tablj2]). The M nn distributions for proton- and deuteron-induced re- 
actions is presented in Figfl^l These distributions are rather wide with a pronounced 
peak in the region of (2-3)m 7r . The fraction of interactions identified with the forma- 
tion of the i?-resonance is slightly above 20 % of all nir collisions. Interactions with 
M n7T ~ 650 MeV resulting in heavy mesons p,u,r] come to about 5 %. 

If the R resonance has been formed, it is assumed to decay only in two photons. 
This scheme can be easily realized by the Monte Carlo method within our transport 
model. The two-photon invariant mass spectra calculated with the inclusion of the 
possible R production are compared with experiment in Fig{TJ3 Indeed, the essential 
part of 77 pairs survives through the strict selected rules and can explain (20 — 30)% 
of the enhancement in the case of dC collisions. For pC collisions the R contribution 
is smaller but in agreement with experimental points. Nevertheless, one should be 
careful in taking too seriously the absolute values of the 77 pair yields estimated from 
the R decay. They are obtained under extreme assumption that all R resonances do 
decay via the two-photon channel, i.e. = T 77 . However, the scalar resonance (like 
the u meson) decays mainly into two pions, Tf> = T n7T , and the electromagnetic decay 
is strongly suppressed. The two-photon decay is dominating (r^ = T 77 ) only if the 
R mass is below the two-pion threshold. So the proposed mechanism allows one to 
consider properly the kinematics of the 77 pairs and the role of acceptance, but it is 



not able to describe the absolute R yield which should be by a few orders of magnitude 
lower than that presented in FigJT5l 




Figure 16: (color online) The calculated 77 invariant mass distribution in pC (left) 
and dC (right) collisions for selected events with iV 7 = 2. Contributions of different 
channels are shown. Symbols near curves describe sources which photons originate 
from. 

Model simulation allows us to disentangle different 7 sources to clarify the produc- 
tion mechanisms. This is illustrated in FigJTBl The double symbol near curves indicates 
the sources where both photons came from. The n° decay (marked as "7r° — 7r°" in 
figures) naturally dominates in both the reactions. In the energy range considered, 
the pion yield rapidly increases and due to that a number of 77 pairs is higher in the 
pC(4.6 GeV) than in the dC(2 AGeV) collisions. The 77 decay ("77 — 77") is seen clearly, 
being spread essentially due to uncertainties in the 7 energy measurement. In the dC 
case the 77 maximum is more pronounced in the total distribution. It is of interest 
that the R resonance decay (" R — > 27") is also visible under PHOTON-2 conditions. 
A number of 77 pairs from R is higher in the pC case, but the ratios for r]/R are 
comparable: t]/R= 34.18 and 24.95 for dC and pC collisions, respectively. It means 
that a possibility to observe the R resonance depends on statistics of measured events. 

One should note that since the low-mass enhancement in the invariant ttti spectra 
showed up clearly at beam energies corresponding to the excitation of A's in the nuclear 
system, the ABC effect was interpreted by a A A excitation [IBJIlT]. In particular, the 
early simplest model for ABC production in pn — > d + X° involves the excitation of 
both nucleons into A-isobars through a one-pion exchange where, after the decay of two 
A's , the final neutron-proton pair sticks together to produce the observed deuteron 
|46j . Though the enhancement observed in the inclusive data for the 7r°7r° channel 
turns out in some cases to be much larger than the predicted in these calculations, 



the AA mechanism is still attractive. More delicate results on the vector and tensor 
analyzing powers in dd — > 4 B.e+X give strong quantitative support to this idea [7] 

The channel marked in Fig {16] as "7— A" corresponds to the case when photons from 
the A decay correlate with any other. Though the two-photon yield for this channel 
in the dC case is close to that from the R decay, the maximum location is shifted to 
higher M 77 by more than 100 MeV. If one chooses both the photons from different A 
isobars, for the PHOTON-2 conditions we have none event from 10 9 of simulated ones. 
One should note that we consider incoherent AA interactions and possible attraction 
in this channel is not taken into account. 

4.3 The 77 -> 3tt° decay 

As was indicated many years ago, a three-pion resonance having the mass M v =550 
MeV and isospin T = maybe a pseudoscalar particle of positive G parity (0~ + ). 
Under this assumption it was shown that the partial rates and widths for this rj decay 
will be consistent with available experimental data providing that the 77 — > 3n channel 
is enhanced by a strong final state interaction [H]. This strong interaction is realized 
by postulating the existence of a particle having the spin and parity + and called a 
a. Then the 3ir decay of the 77 meson would proceed in two distinct steps: 77 — > a + 7r°, 
a — > 2tt° or (7r + 7r - ) where the first step is an electromagnetic decay, while the second 
occurs through the strong interaction. The fit to experimental data gives the mass of 
a cr particle about 370 MeV and a full width of about 50 MeV [H]. These parameters 
coincide with those extracted from the direct analysis of pion spectra from the 3ir 
decay of the 77 meson [19J [50]. The enhancement of the 3ir channel was argued by 
the presence of a strong two-pion interaction which resembles the ABC effect. It is of 
interest that the presence of such a pion-pion resonance improves also the calculation 
of the K L — K s mass difference [51] . 

To check this mechanism (see Fig JTTj) we simulated two channels of the 77 decay: 
the direct decay into two photons 77 — > 27 and t] —> 37r° which then decay into photons. 
The last channel was calculated under two assumptions. The first, all the three pions 
decay independently, 7r° — > 77, creating 6 photons. The second version assumes that, 
as discussed above, two pions may interact forming a which decay into 2 pions, so 
77 — > a + 7T° — ► 67. As it is seen (bottom panels in Fig JT7|) . the interaction in the inr 
channel results in some enhancement in the ttti invariant mass spectra as compared the 
noninteracting case. The 77 invariant mass spectra exhibit a spread maximum near 
the pion mass and they are practically identical in both the cases (dashed lines). If the 
PHOTON-2 selection conditions are implemented, a clear signal (solid lines in Fig JTTj) 
at M 77 ~ 350 MeV appears but its intensity is by about three orders of magnitude 
lower than that for the 77 meson. It is of interest that the absolute values and shape 
of the M 77 spectrum are again very similar in both the versions. So in the 77 — > 3tt° 
decay it is hardly ever possible to disentangle the cases with and without the two-pion 
interaction by the detection of decay photons. 

4.4 Dibaryon mechanism 

Recently, a resonance-like structure has been found by the CELSIUS-WASA Collab- 
oration in the two-photon invariant mass spectrum near M 77 ~ 2m n of the exclusive 
reaction pp — > 727)77 at 1.2 and 1.36 GeV [52]. This observation was interpreted as a 
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Figure 17: (color online) The 77 (top) and ttit (bottom) invariant mass distributions 
for the 77 decay through the 77 — > 77 and 77 — > 37r° channels from pp collisions. The left 
panels correspond to the 77 — > 3tt° — > 67 mechanism, results in the right panels include 
interaction of two pions via the a meson, 77 — > a + it — * 67. The dashed lines are 
calculated for the full Air acceptance, the solid lines take into account the PHOTON-2 
kinematical conditions. 

channel pion loops which are generated by the pp collision process and decay into the 
77 channel. Interference with the underlying double bremsstrahlung background can 
give a reasonable account of data [52]. 

In Ref. [53], some arguments were given that such an interpretation is at least ques- 
tionable and an alternative explanation was proposed where a possible origin of the 
structure is based on the dibaryon mechanism of the two-photon emission [51]. The 
proposed mechanism NN d\ ^ NN'-fy proceeds trough a sequential emission of 
two photons, one of which is caused by production of the decoupled baryon resonance 
d* and the other is its subsequent decay. The pp — > pp'yy transition is treated in [53] 
within the assumption that at a large distance the NN- decoupled six quark d\ state 
is a bound pA(1232) state with the spin-parity J p = 0~ and isospin 1 = 2 [55J. The 
matrix elements were estimated phenomenologically and effects of the final state inter- 
actions between decay protons were included. This model reproduces reasonably well 
the experimentally observed M 77 spectrum of the pp — > pp'yy reaction in the vicinity 
of the resonance structure [53] . 

We would like to check whether this dibaryon mechanism may be responsible for the 
peak observed in the M 77 distribution of the dC collisions at T = 2 AGeV. The two- 
step scheme iViV — > ri^7 — > ./V./V77, where the dibaryon mass is = mN + m/\ = 2.182 
GeV, can be easily simulated and included into our transport model. The only unknown 
quantity is the cross section of this process. We estimated it by means of the linear 




Figure 18: (color online) The 77 invariant mass distribution (left) and energy spectra of 
photons (right) calculated for dC collisions with inclusion of the dibaryon mechanism. 
Contributions of different channels are shown similarly to FigJTBl All curves beside the 
dibaryon one are given for the PHOTON-2 selected events. For the dibaryon channel 
(marked as "dibar") the total two photon yield is presented. 

extrapolation of the two available points measured at 1.2 and 1.36 GeV [52] till the 
energy of about 2 GeV. 

In FigiTBl such model calculation results with inclusion of the dibaryon channel are 
presented. If the beam energy of the pp collision is fixed by 2 GeV, the photon energy 
at the first step pp — > d\^ will be a line, E y i = (s — M%)/(2y/s) = 640 MeV in the 
cm. system. Here s is the total pp colliding energy squared. In the lab. system the 
maximal photon energy reaches about 1.5 GeV. However, the photon from the second 
stage will be soft, being defined by the baryon mass. Nevertheless, the maximum 
position of appropriate two photon distribution moves also with the energy increase. 
The total photon spectra from dibaryon and other sources from dC interactions are 
shown in the right panel of Fig (18] and it has a two-bump structure which mirrors 
the two-step production mechanism. One should note that in contrast with all other 
results in this figure the distributions of the dibaryon channel are obtained for the full 
4n acceptance without any limitation on energy. If the PHOTON-2 selection conditions 
are implemented to the dibaryon channel, the low energy photons are cut and we get 
no two-photon pair from the dibaryon among 10 9 simulated collisions. Therefore, this 
dibaryon mechanism cannot explain the observed anomaly. 

One should note that in a certain sense the model considered is a conventional 
dibaryon model where no nonhadronic degrees of freedom are involved. Generally, 
these results may differ from those obtained within nonconventional dibaryon models 
like ESI. 



5 Estimate of the r\ and R production cross sections 
and resonance widths 



The summed number of pC- and c/C-interactions in the experiment amounts to ~ 3-10 11 
and ~ 2 ■ 10 12 , respectively. The inelastic cross sections of the observed pC and dC 
reactions are u inc i(pC) =411 mb and a ine i(dC) = 426 mb [57], respectively. 

The cross section for the r] production in dC collisions (similarly to pC interactions) 
is defined as follows: 

/yexp ?ymod 

a(dC ^ V + X) = a incl (dC) ■ Ar f)l . (8) 

^dC-intcr N™ d / K opt 

Here the first two factors are the reaction cross section and the measured mean 
multiplicity of t] mesons. A number of true dC interactions resulting in the rj production 
is given as 

ArfC— inter -^-empty ' -^beam— absorb " -^d i (9) 

where the total number of beam particles passing through the target = (l-j-2.5) -10 12 
is corrected on possible interactions outside the target i^empty ~ 0.995, to be estimated 
by a special experiment with an empty target, as discussed in Sect. 12.21 and the beam 
absorption .^beam-absorb = 0.5 ± 0.2. The last factor reflects a particularity of internal 
beam experiments where the direct monitoring is impossible and part of the initial 
beam does not interact during the working circle. 

The third factor in eq.© is the ratio of the total number of simulated 77 mesons to 
a number of 77's decaying into two photons under PHOTON-2 experimental conditions. 
The coefficient K opt = 8.93 takes into account the rotation of the modelled events in 
the 4> plane to find other possible 77 pairs in the given event satisfying the selection 
conditions. Assuming that photons in the event are distributed homogeneously in 0, 
this trick allows one to increase effectively statistics of the selected events. 

So for the 77 production in dC collisions we have 

a(dC ^ 77 + X) = 1.31 ± 0.11±5;|| mb . 

In the case of pC collisions we get N p = (1.5 -j- 5) • 10 11 , i^beam- absorb = 0.8 ± 0.2, 
coefficients K empt y and K opt are practically the same as for <iC-collisions, so for the 
cross section we have 



a(pC ^ rj + X) = 3.2 ± 0.2j£jj mb . 

Large systematic errors of cross sections are coming mainly from a problem of moni- 
toring the intensity of the internal beam. 

The obtained values for cross sections are compared in Fig{19] with calculated ex- 
citation functions for the rj meson production in the pC reaction. Note that the DCM 
model agrees in the absolute scale with the experimental differential cross section for 
this reaction in the energy range of T p = 0.8 -j- 2.0 GeV (see Fig JTT]) . The measured 
c{pC — > rj + X) is below theoretical predictions by a factor of about 2. The scaled 
pp excitation essentially differs from nuclear one in the near-threshold region due to 
Fermi motion. The needed scale factor A = 12 is a little bit higher than naively ex- 
pected A 2 / 3 because these data correspond to the pp — > ppr] channel only while nuclear 




Figure 19: (color online) Excitation function for 77 production in pC collisions calculated 
by Cassing [SB] (dotted curve) and within our DCM model (solid line). Our measured 
points for dC (with a factor of 0.5) and pC collisions are plotted by the star and 
square, respectively. Circles show the cross section for elementary collisions pp — > ppr\ 
multiplied by the factor of 12. 



excitation functions include all channels, pp — ► rj + X. Unfortunately, there is no other 
experimental points in this energy range. 

If the cross section for the 77 production is known, the cross section for the R- 
resonance may be estimated as follows: 

N ex P(R -> 77) leu 

a(dC^R^ 11 )=a(dC^, + X) • Br(„ - 77) • N J {v ^y* = 
= (0.075 ± 0.018) • a{dC -> 77 + X) = 98 ± 24^ fib , (10) 

where the branching ratio is Br(r] -» 77) = 0.38, e R = N mod (R -» ^)/N™ d (R) 
and e v = N mod (r] — > 77)/-?V t 1 ot od ( ? 7) are ^ ne detection and selection efficiency The mea- 
sured average reduced multiplicities are compared in Figj20]with available systematics 
for meson production near the threshold energies. This simple energy scaling system- 
atics for the subthreshold and near threshold particle production was proposed in [BU] • 
A number of participant nucleons is estimated from the geometrical consideration as 



A n A 2 ^ + A T A 2 J 3 
(Ap + Aj 



P I A 1 / 3 J /lV3\2 



The bombarding kinetic energy is corrected for the Coulomb barrier V c . This system- 
atics is valid also for K and p mesons [591 EH] . As is seen, the measured r\ production 
points are somewhat below the general trend. Partially, it may be caused also by that 
all other experimental points correspond to heavier nuclear systems and eq. ( fTTI) is not 
very justified for our reactions. 
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Figure 20: (color online) Meson production probability per participant nucleon as a 
function of the bombarding energy per nucleon normalized to the respective meson 
production threshold T t h r in free nucleon-nucleon collisions. Experimental points of 
the TAPS collaboration are taken from the review [59]. Our points for dC and pC 
collisions are plotted by the star and filled square. 



As to the true internal width w of the observed resonances, they are defined by the 
measured width w me 



and also specified by the spectrometer resolution w sp 

1 1/2 



W 



w 



spy 



(12) 



The spectrometer resolution depends on the M 77 range. For the R and r] invariant 
mass range we have, respectively, 



2 w; S p (340 < M 77 < 360MeV) 
2w sp (540 < M 77 < 560MeV) 



52.6 MeV, 
68.6 MeV. 



So, according to eq.l jlzj) . the intrinsic widths of detected resonances are 



(13) 



2w(rj 
2w(R 



77) 
77) 



0. 

63.7 ± 17.8 



(14) 



As is expected, the width of the ?7-meson practically equals 0, whereas it essentially 
differs from zero for the observed resonance. The value of 2w in the Gaussian distri- 
bution (J3]) practically coincides with the width T in the Breit-Wigner function; thus, 
the intrinsic width of the observed resonance structure is about 64 ± 18 MeV. 



6 Concluding remarks 

Thus, based on a thorough analysis of experimental data measured at the JINR Nu- 
clotron and statistics of 2339 ± 340 events of 1.5 • 10 6 triggered interactions of a record 
total number 2-10 12 of dC- interactions there was observed a resonance-like enhancement 



at the mass M 77 = 360±7±9 MeV, with the width r = 64 ±18 MeV. The production 
cross section a 77 ~ 98 fib is estimated preliminary in the invariant mass spectrum 
of two photons produced in dC- interactions at momentum of incident deuterons 2.75 
GeV/c per nucleon. A structure like this is not observed in the M 77 spectrum from pC 
(5.5 GeV/c) interactions while the r\ meson is clearly seen in both the cases. These re- 
sults, obtained by means of the mixing event background, are confirmed by the wavelet 
analysis. 

Due to the use of internal beams of the JINR Nuclotron, totally more than 10 6 
triggered events at high discriminator thresholds were recorded for every reaction dur- 
ing these experiments. Smallness of the signal-to-background ratio in the R-resonance 
mass range needs very high statistics for observation. This demand was not satisfied 
in previous experiments and therefore no structure has been observed in this invariant 
mass range. As was noted above, the only most statistically meaningful measurement 
of the invariant mass spectra was made by the TAPS-collaboration but to resolve the 
resonance structure discussed, a number of registered events in the TAPS experiment 
should be increased by an order of magnitude. 

To certain extent this enhancement at M 77 ~ (2 — 3)m. K is similar to the puzzling 
ABC effect observed for two-pion pairs from nucleon-nucleon and lightest nuclei col- 
lisions at the near threshold energy. In the given work we see that it exists in the 
77 channel and measurements are extended to a heavier system. It means that this 
resonance-like structure is a quite stable object which even survives in the nuclear 
surrounding. 

To understand the origin of the observed structure, several dynamic mechanisms 
were attempted: production of the hypothetic R resonance in inr interactions during 
the evolution of the nuclear collision, formation of the R resonance with participation 
of photons from the A decay, the 7r°7r° interaction effect in the 37r° channel of the 77 
decay, and a particular decoupled dibaryon mechanism. Unfortunately, none of these 
mechanisms is able to explain the measured value of the resonance-like enhancement, 
though they contribute to the invariant mass in the region of interest. 

The carbon target is really the heaviest one used in experiments where ABC-like 
structure has been observed. In contrast with all other experiments considered here, 
one may expect some manifestation of in-medium effects. The prominent feature of the 
77 meson is that the ?7-nucleon system couples dominantly to the A r *(1535)(S' 11 ) reso- 
nance at the threshold energies. Hence, due to the 77 coupling to A r *(1535)-nucleon-hole 
modes, one could expect the eta meson nuclear dynamics to be sensitive to modifica- 
tion of nucleon and N* properties in medium. As was shown in [BT], the 77 spectral 
function at normal nuclear density has a second maximum near m v ~ 400 MeV which 
may be associated to a partial chiral symmetry restoration. Its two-photon decay 
inside a nucleus might give a rise to a maximum close to the measured value of R. 
Unfortunately, we cannot perform transport calculations with taking into account the 
in-medium modification of hadron properties. 

As was noted in Introduction, the recent data of the WASA-CELSIUS Collaboration 
provide a strong support to the idea of a nontrivial dibaryon state [171 ITS] . An attractive 
candidate for its realization may be a model of the intermediate a-dressed dibaryon [56] • 
In this model the short-range NN interaction, described by the standard t-channel a 
exchange between two nucleons, is replaced with the respective s-channel a exchange 
associated with the intermediate dibaryon production treated as a <r-dressed six-quark 
bag. The strong scalar a-field arises around the symmetric 6q bag, because the change 



in the symmetry of six-quark state in the transition from the NN channel to the 
intermediate dressed-bag state. Due to a strong attraction of the a meson to quarks, 
this intensive a field squeezes the bag and increases its density. The contribution of the 
s channel mechanism would be generally much larger than the conventional t-channel 
one due to a resonance-like enhancement. The high quark density in the symmetric 6q 
state enhances meson field fluctuations around the multiquark bag and thereby partially 
restores the chiral symmetry. Therefore, the mass of a meson gets much lower and has 
been estimated to be the value m a ~ 350 -r- 380 MeV. In its turn, it should enhance the 
near-threshold pion and double-pion production [56], [62]. In addition, a large variety 
of nuclear data, in particular properties of short- and intermediate-range of NN and 
3N potentials, has been explained within this model; however, still there is no direct 
quantitative calculations of the ABC-like effects. 

From the experimental side it is highly desirable to determine more accurately 
the mass, width, and cross section of the observed resonance structure by enlarging 
the acceptance. To verify the above conclusions new experiments are required to be 
carried out under conditions appropriate for detection of pairs of two photons within 
the invariant mass interval of 300-400 MeV. In this respect experiments on proton and 
carbon targets with proton and deuteron beams at the same energy per nucleon would 
be very useful. Some scanning in the beam energy will clarify the possible resonance 
structure of this effect. By varying the opening angle of the PHOTON-2 spectrometer it 
is possible to get information about momentum spectra of the produced resonance- like 
structure which could be a delicate test of the R production mechanism. 
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Appendix A. Continuous wavelets with vanishing moments Here we would 
like to elucidate some details of the wavelet analysis. 

The family of continuous wavelets with vanishing moments (VMW) is presented 
here by Gaussian wavelets (GW) which are generated by derivatives of the Gaussian 
function ([7]). For canonical Gaussian with xq = 0; o = 1 and A = 1 one obtains 

Jn 

G n (a, b) = G n (x) = (-l)«+^ e -* 2 / 2 , (15) 

ax n 

where n > is the order of the g n {x) wavelet. The normalizing coefficients of these 
wavelets C 9n are 2n(n — 1)!. 



The most known in the GW family is the second order GW 



G 2 (x) = (l-x 2 )e-^, 



which is also known as "the Mexican hat" [3Tj . 

We use here also GW of higher orders, in particular, 



„.2 

„2 „4 



G 4 (x) = (6x z - x 4 - 3)e"- (16) 
G 6 (z) = -(x 6 - 15x 4 + 45x 2 - 15)e - T (17) 



„2 



G 8 (x) = -(x 8 - 28x b + 210x 4 - 420x 2 + 105)e"- (18) 

It is a remarkable fact that the wavelet transformation of Gaussian (J7J looks like the 
corresponding wavelet. Therefore, the general expression for the n-th wavelet coefficient 
has the following form (see derivations in |30j): 

W gn (a,b)g = g n — ° , (19) 

y/{n- 1)! s n+1 \ s J 

where we denote s = y/ a 2 + a 2 . 

The parameters a and b of continuous wavelets in eq. (jlj) are changing continuously, 
which leads to the redundant representation of the data. In some cases, the above 
mentioned GW properties are quite useful, in particular, this redundant representation 
facilitates careful spectrum manipulations. The price of this redundancy consists in 
slow speed of calculations. Besides, all signals to be analyzed have in practice a discrete 
structure. 

It is noteworthy that GW should be used with some care since they are nonorthog- 
onal, which may disturb amplitudes of the filtered signal after their inverse transform. 
In this respect the discrete wavelet transform looks more preferable for many applica- 
tions of computing calculations with real data [32] and deserves special consideration. 
As was noted above, in our particular case of the continuous WVM one can identify 
resonances without the inverse transformation. 

Let us demonstrate this scheme by the £2(0, b) wavelet example. 

According to eq. f|T9|) . one can obtain the maximum (absolute) value of G2 for a 
Gaussian (J7|) at the shift point b = Xq as 

Aaa 5/2 

max W G2 {a,x )g = — rr^. (20) 

In the wavelet domain of (^2(0, Xq) this dependence looks like a simple curve with one 
maximum. To find it, one has to solve the equation 

^- max W G2 (a, x )g = 
da b 

The corresponding calculations give the position of maximum at the scale 
a-max = V^c Since the maximum location in the G2 domain is stable when the signal 
is contaminated by some noise (see [30]), one can use the obtained point xo,a max to 
start the fit. Although the maximum of a real contaminated signal is inevitably blurred 
over some area in the wavelet space due to various distortions, it can nevertheless be 
used as a good starting point for iterations minimizing a fitting functional. 
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